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Abstract 

We show how to use the lattice Green function to calculate capacitances in two dimensions with bound- 
ary conditions at infinity. It is shown how to calculate coefficients of capacitance and induction from the 
lattice Green function. A general analysis of two arbitrary conductors is carried out. It is shown how the 
calculations can be simplified in the case of identical conductors when certain symmetry conditions are met. 
Example calculations for a parallel and coplanar stripline are shown. The use of the two conductor formulas 
for the case of three or more conductors is discussed. 



*Electronic address: rstefan@exstom.coml URL: |http : / /www, exstrom. com/stef an/stef an .html | 



1 



I. INTRODUCTION 



In a two dimensional homogeneous and unbounded space, the Poisson equation and its solution 
are given by 

V^(?) = -P© (1) 

* (?) = 2k/ 1 "(^) p(? ' ,< "' <2) 

When the charge density is known then eq. reduces the problem of finding the potential at 
any point in space to a simple integration. A more interesting and difficult problem occurs when 
the charge density is not known. An example of such a probelm is the case of two conductors 
each of which is held at some constant potential. Eq. |2]is then an integral equation for the charge 
density on the conductors and it is difficult to solve analytically in all but a few simple cases. Many 
approximation schemes have been devised to deal with this problem and in most cases the solution 
must be found numerically. One commonly used approach is to discretize eq. El thereby turning it 
into a matrix equation, which can then be solved using standard techniques. The disadvantage of 
this approach is that it is an approximation of the continuous as well as the discrete version of the 
problem, i.e. it is an approximation in both continuous and discrete space. 

An alternative approach is to formulate the entire problem in a discrete space. Physically the 
model is an infinite square lattice with capacitors connecting the nodes. Sets of adjacent nodes, that 
correspond to discretized versions of conductors, are held at constant potential and the problem 
is to find the charges at those nodes. Mathematically this approach is equivalent to replacing the 
Laplacian in eq. [TJwith its finite difference approximation. The discrete space version of eq. [T]is 
then 

Y, L nmWm) = ~— — (3) 

m ^ 

where r n is a lattice vector of the form: r n = n\a\ + njai, «j =integer, a ; = axi and a =lattice 
spacing. \{r n ) is the linear charge density at node r n . L nm is a matrix element of the lattice 
Laplacian and is defined as 
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Lnm = -48(r„, r m ) + £ [8(r„ + a,-, r m ) + 8(r„ - fl/, r m )} (4) 
i=l 

To simplify the notation we will write eq. |3]in the following form 

Y, L nm§(r m ) = -qif n ) (5) 



where L nm now has units of capacitance and q(r n ) is a charge at node f n . The solution of this 
equation for §(r n ) is 

ty{r n ) =Y i G mn q(r m ) (6) 

m 

where G nm is a matrix element of the Green function with units of elastance. The Green function 
in this problem is defined as, G = —L , where L is the lattice Laplacian operator. Eq. |6]is the 
discrete space version of eq. |2j 

G nm is a function of \f n — r m \ only, so a more convenient notation for eq. |6]is 

§(m,n 2 )= G{p u p 2 )q(m u m 2 ) (7) 

mi,ni2 

where p\ = \n\ — m\\ and p 2 = \n 2 — m 2 \. The problem with this equation however, is that 
G(p\,p 2 ) is infinite for all values of p\ and p 2 . This same problem occurs in the continuous 
case, eq. \2\ when the reference point for the potential, ro, is taken to infinity. The way around this 
problem is to use g(pi,p 2 ) = G(0,0) — G(pi,p 2 ) which is finite for all values of p\ and p 2 
Eq. |7]then becomes 

§(n h n 2 ) = - £ g(p u p 2 )q(m u m 2 ) (8) 

m u m 2 

This equation will solve the original problem as long as all the charges sum to zero. We will now 
show how to use this formalism to calculate the capacitance between conductors in two dimen- 
sions. 



II. GENERAL ANALYSIS OF TWO CONDUCTORS 



Consider the case of two conductors discretized so that there are n\ charges on conductor 1 and 
n 2 charges on conductor 2. The conductors are held at constant potentials (|)i and § 2 - The equation 
for this system is written as follows 

/ Gu Gn W qi \ / <tiii \ (9) 
yG 2 i G 2 2 J \q2 J J 

q ; is an n ; dimensional vector of the charges on conductor i. e, is an n, dimensional vector with all 
elements equal to 1. Gu is an matrix that gives the contribution of the charges on conductor 
i to the potential of the conductor. G// is an nixrij matrix that gives the contribution of the charges 
on conductor j to the potential of conductor i. The elements of G (J depend only on the absolute 
separation of two charges therefore it will always be true that G,- ; = Gj). 
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As a simple example, take two conductors discretized such that there are two charges on both 
conductors. Let the charges on conductor 1 be at (0,0) and (0, 1) and the charges on conductor 2 
beat (2,0) and (3,0). The G (J submatrices in eq. [5] will then be 

I g(l,0) 

G n =G 2 2= | (10) 

\s(i,o) o 
\«(2,1) «(3,i; 

In general, to calculate the capacitance between two conductors the matrix in eq. [9]needs to be 
inverted in order to find the charges qi, and q2. The matrix is easily inverted if it is first factored 
into a product of an upper and lower triangular matrix. This can be done in two ways 

I o\ / Gn G12 

G21G11 I / \ G22 — G21&J/G12 



or 





(12) 



I Gi2G 2 2 \ / Gn — Gi2G 22 1 G2i 



(13) 



I / \ G21 G22 

Inverting these upper and lower triangular matrices then leads to two sets of equations for the 
submatrices of the inverse. From eq. [12] we get 

(G- 1 )!! = Gr 1 1 + Gri 1 Gi2(G 2 2-G 2 iGr 1 1 Gi2)- 1 G2iGr 1 1 (14) 

(G l )n = — G n 1 Gi2(G22 — G21GJ/G12) 1 

(G _1 )2i = — (G22 — G2iG^Gn) 1 G2iG 11 1 
(G )22 = (G22 — G2iG u 1 Gi2) 1 

and from eq. [13] we get 

(G -1 )n = (Gn-GijG^Gji)- 1 (15) 

(G _1 )i2 = — (Gn — Gi2G 22 1 G2i) G12G22 

(G _1 )2i = — G 22 G2i(Gn — Gi2G 22 G2i) 1 

(G ! )22 = G 22 +G 22 G2i(Gn — Gi2G 22 G2i) G^G^ 

In terms of these submatrices the charges on the two conductors are given by 

qi = -^i(G- 1 )nei-^2(G- 1 )i2e 2 (16) 
q 2 = -4>i(G _1 ) 2 iei -4> 2 (G _1 ) 2 2e2 



The total charge on conductor i is Qi = efq ; , so if the first equation is multiplied by ef and the 
second equation by men we g et a set °f equations relating the total charges on each conductor 
to their potentials. 

Q\ = c n §i+ci2§2 (17) 

The coefficients cu are known as coefficients of capacitance and c\i is known as a coefficient of 
induction. The coefficient c, 7 is equal to the negative of the sum of all the elements in the matrix 

c ij = -eJ(G- l ) ij ej (18) 

From eq. [14] and [15] it is clear that (G _1 ) i; - = (G -1 )^- and therefore Cij = cji. Since the potential 
and the charge of a conductor will have the same sign, we have the condition cu > 0. Also since 
the charge induced by a conductor will have a sign opposite to the potential of the conductor, we 
have cu < 0. 

In addition to eq. [17] we have the requirement that the sum of all the charges must equal zero, 
Qi + Qi = 0- This means that we can set Qi = Q and Q2 = —Q in eq. [T7]and solve for the 
potentials. 

♦. - ( C22+C12 2 )g < 19 » 

\C\\C22-C\ 1 ) 
. ( C U +C12 \ 

<t>2 = - Y Q 

\CUC22-C\ 2 J 

The ratio of the two potentials is then given by 



\ , - (2°) 

<P2 \ c ll+ c 12, 

The capacitance between the two conductors can then be expressed in terms of the Cjj coefficients 
as follows. 

c= Q _ cnc22-c 2 12 

§l-§2 Cn+C22+2ci2 

Note that the denominator, c\ \ + C22 + 2c 12 is equal to the negative of the sum of all the elements 
of the inverse of the Green function matrix in eq. [5] and that the c ;/ 's are functions only of the 
geometry of the two conductors. Eq. ET]can be used to calculate the capacitance between two 
conductors of arbitrary size, shape and orientation once the ctj coefficients have been calculated. 
Each conductor may consist of more than one disconnected piece as long as each piece is held 



at the same potential. We now look at the special case of two identical conductors, i.e. both 
conductors have the same size and shape. 



HI. TWO IDENTICAL CONDUCTORS 

For two conductors of the same size and shape the calculation of capacitance and charge distri- 
bution can be simplified in those cases where the coefficients of capacitance are equal. We begin 
then by examining what conditions are required to get en = 022- 

Since cu is equal to the sum of the elements of the (G _1 ) ;; - matrix, the two c« will automatically 
be equal if the two (G _1 )/ ; matrices are equal. For identical conductors it will always be true that 
Gn = G22 so that the expressions for (G _1 ),- ; in equations I"R1 and [T31become 

(G- 1 )!! = (Gn-GizGr/Gzir 1 (22) 
(G )22 = (Gn — G2iG^ 1 G12) 1 

These expressions are only equal if G12 = G21 = &\ 2 G12 must be symmetric. 

G12 is a function only of the distance between charges on the two conductors therefore it can be 
made symmetric if it is possible to take one conductor to the other by a reflection about a horizontal 
plane, a vertical plane, both a horizontal and vertical plance, or a diagonal plane. Figure 1 shows 
the case of two conductors that can be made congruent through reflection about the vertical plane 
ab followed by reflection about the horizontal plane cel. The point 1 is taken to the point 1' so 
that these charges will be equal and all charges to the right of point 1 will equal the corresponding 
charges to the left of point 1'. 

With G12 symmetric we have c\\ = C22 and the equations for the general case simplify. From 
eq. |2Hl we get §\ = — §2 = § and the potential is related to the charge as follows 

4> = Q (23) 
cil -C12 



The capacitance is then given by 



C= C -l^ (24) 



For this symmetric case, each charge on one conductor will have a corresponding equal and oppo- 
site charge on the other conductor so that in eq. |9]we get qi = — q2 = q and the equation reduces 
to 

(Gi2-Gii)q = ^e (25) 



The total charge on one conductor is then related to the potential as follows 



2 = ^e r (Gi2-Gn)- 1 e 



(26) 



Comparing this equation with eq. 1231 gives 



c\\ -C12 = e r (Gi2-Gn) 'e 



(27) 



The capacitance can therefore be calculated from the sum of all the elements of the inverse of the 
matrix G12 -Gn. 

As examples, we will now consider the case of the parallel stripline and the coplanar stripline. 
A parallel stripline and the potential surrounding it is shown in Fig. 2. The upper and lower plates 
are at potentials +1 and -1 respectively. A plot of the capacitance as a function of the ratio of 
the width of the plates to their separation is shown in fig. 3. For each ratio the resolution was 
increased (lattice constant decreased) until the capacitance value appeared to converge. Fig. 4 
shows the convergence for the ratios 1, 5 and 10 as a function of the number of charges in the 
plates. With r equal to the ratio of the width to the separation, the following equation fits the plot 
of the capacitance in fig. 3 with a correlation coefficient of 0.999998. 



This equation agrees well with previous work by Wheeler [3]. The coplanar stripline and the 
potential surrounding it is shown in fig. 5. The left and right plates are at potentials -1 and +1 
respectively, with a plate separation of 1/5 the width of a plate. A plot of the capacitance as a 
function of the ratio of width to separation of the plates is shown in fig. 6. 

IV. CONCLUSION 

The formalism discussed above can easily be extended to the case of three or more conductors. 
The Green function matrix inverse formulas in eq. [T3]and[T3]can be applied iteratively in this case. 
Start with any two conductors, generate G\\, G22, G\2, and then use one of the sets of equations 
to find the inverse. This inverse then becomes G u ' in eq. fl4lwhen the next conductor is included, 
for which we have a new G22 and G\2, with G\2 connecting the new conductor with the previous 
two. This makes it possible to look at the effect of changes in the placement of a new conductor 
(or a single charge) with respect to a set of conductors without major recalculation efforts. 



C = r+ 1.12863 + 0.202069 ln(r- 0.0964012) 



(28) 
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Another theorem that may be useful in the case of three or more conductors is Green's Recip- 
rocation Theorem. In its simplest form the theorem says that if we have a set of N conductors at 
potentials and charges Qi (i = 1, . . . AO and then take those same conductors at new potentials 
$ and charges Q\ then the following equation holds 

£e^; = £e>/ (29) 

1=1 7=1 

This theorem is easily proven by writing out equations similar to eq. [T7]for the Qi and Q' t . The 
equations for Qi are multiplied by and the equations for Q^ are multiplied by fy. Summing the 
two sets of equations then gives eq. 1291 

The capacitance calculations we have discussed can also be applied to finding the energy and 
forces on conductors and to calculating transmission line impedances. The energy of a set of N 

conductors with charges Qi and potentials (|), is |Q] 

i N i N N 

W = -Y t Q i to = -Y t Y.c u M J (30) 

!=1 i=l 7=1 

To calculate the force on a conductor in a given direction, with all conductors held at constant 
potential, we displace the conductor in that direction and then calculate the new c'^s and the new 
energy. The force is then the change in energy divided by the lattice constant. 

The impedance of the transmission line formed by two conductors is given by the equation 

Z=£ (31) 

where T] is the characteristic impedance of the medium, given by T| = ^/fi/e, and C is the ca- 
pacitance given by eq. |2T] Note that C is treated here as a pure number i.e. it is not scaled by 
e. 

Formulating electrostatics problems and their solution entirely in discrete space has many ad- 
vantages. This was first recognized in a well known paper by Courant et al |Q] in which they 
examined the solution of elliptic partial differential equations and their corresponding difference 
equations. They were able to show that the difference equation solution does converge to the 
solution of the differential equation and that some questions, such as the existence of solutions, 
are more easily answered by looking at the difference equation. It appears that in many cases 
the theorems and methods developed for continuous space problems can be translated over into 
discrete space. It also seems to us that the possibility exists for discovering new theorems or gen- 
eralizations of existing theorems, such as Thompson-Lampard [6], by approaching problems from 
a discrete space point of view. Opportunities for more research in this area certainly exist. 
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Figure 1 : Symmetry example. 




Figure 2: Parallel stripline potential: 100 charges/plate, separation is 1/2 plate width. 
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Parallel Stripline Capacitance 
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Figure 3: Parallel stripline capacitance: F/m in units of epsilon vs. width/separation. 1000 charges per 
plate. 
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Figure 4: Parallel plate convergence: F/m in units of epsilon vs. charges per plate 
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Figure 5: Coplanar stripline potential: 100 charges/plate, separation is 1/5 plate width. 
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Figure 6: Coplanar stripline capacitance: F/m in units of epsilon vs. separation in units of lattice constants. 
300 charges per plate. 
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